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ABSTRACT 

We present a power spectral analysis of a 100 ksec XMM-Newton observation of 
the narrow line Seyfert 1 galaxy Ark 564. When combined with earlier RXTE and 
ASCA observations, these data produce a power spectrum covering seven decades of 
frequency which is well described by a power law with two very clear breaks. This 
shape is unlike the power spectra of almost all other AGN observed so far, which 
have only one detected break, and resemble Galactic binary systems in a soft state. 
The power spectrum can also be well described by the sum of two Lorentzian-shaped 
components, the one at higher frequencies having a hard spectrum, similar to those 
seen in Galactic binary systems. Previously we have demonstrated that the lag of the 
hard band variations relative to the soft band in Ark 564 is dependent on variability 
time-scale, as seen in Galactic binary sources. Here we show that the time-scale 
dependence of the lags can be described well using the same two-Lorentzian model 
which describes the power spectrum, assuming that each Lorentzian component has 
a distinct time lag. Thus all X-ray timing evidence points strongly to two discrete, 
localised, regions as the origin of most of the variability. Similar behaviour is seen 
in Galactic X-ray binary systems in most states other than the soft state, i.e. in 
the low-hard and intermediate/ very high states. Given the very high accretion rate 
of Ark 564 the closest analogy is with the very high (intermediate) state rather 
than the low-hard state. We therefore strengthen the comparison between AGN and 
Galactic binary sources beyond previous studies by extending it to the previously 
poorly studied very high accretion rate regime. 

Key words: Galaxies: active, powerspectra: Lorentzian components, accretion: 
discs, spectral timing: phase lags 



1 INTRODUCTION 

Understanding the relationship between AGN, which are 
powered by accretion onto supermassive black holes, and 
the stellar mass Galactic black hole binary systems (GBHs) 
is currently one of the major research areas in high energy 
astrophysics. If we can understand the relationship, then we 
can predict how AGN should behave on cosmological time- 
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scales by observing the brighter, and much faster varying, 
GBHs. 

The comparison is complicated by the fact that (at the 
risk of some oversimplification), GBHs can exist in a num- 
ber of states, defined originally in terms of their medium 
energy (2-10 keV) X-ray properties, particularly: (i) the 
hard spectrum, low-flux ('hard') state, (ii) the soft spec- 
trum, high flux ('soft') state and (iii) the 'very high state' 
(VHS - sometimes also described as the high-intermediate 
state) where the flux is very high and the spectrum very 
soft. In general the accretion rate increases as we go from 
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the ha rd to soft to VHS states. See lRemillard fc McCUntockl 
(|2006l ) for a comprehensive description of GBH states and 
behaviour. 

The various states have quite distinct X-ray timing 
properties which may, in fact, provide a more fundamen- 
tal state discrimator than the spectral properties. Timing 
is usually quantified in terms of the power spectral densi- 
ties (PSDs) of the X-ray lightcurves. Where significant vari- 
ability is observed, 'soft' state PSDs are characterised by a 
power law of slope of -1 at low frequencies (i.e. the power, 
P{v) oc with a ~ 1), breaking to a slope a > 2 above a 
characteristic frequency vb (or time-scale Ts). Hard state 
and VHS state GBH PSDs, which are similar to each other, 
are more complex. They are best described by the combi- 
nation of two or more Lorentzian-shaped components (e.g. 
[Nowak 2000), although at low signal to noise they can be ap- 
proximated by a power law with a second, lower frequency 
break, approximately 1.5 to 2 decades lower in frequency 
than the high frequency break. 

So far, almost all well observed AGN, which are X- 
ray bright and have moderately high accretion rates, have 
soft state PSDs, as we would expect from comparison with 
GBHs of similar accretion rate (me asured relative to the 
Eddi n gton accretion rate , rhg) fe.g. iM'^Hardv et al.ll2004l . 
I2OO5I : luttlev fc M'^Hardvlliooa 'l. However the comparison 
for the hard and VHS states is largely unknown. Until very 
recently, N GC3783 had been sugge sted to have two breaks 
in its PSD (jMarkowitz et alll2003l ) and so be a hard state 
AGN, although its accretion rate is similar to that of the 
well observed soft state AGN. However r ecent work, based 
on improved data l|Summons et al.ll2007l ). has shown that, 
in fact, NGC3783 is another soft state AGN. If soft states 
are linked to the accretion disc extending in to very small 
radii, then perhaps their much cooler accretion discs com- 
pared to GBHs may allow the optically thick disc to survive 
without evaporation to smaller radii. Thus maybe all AGN 
are soft state objects and the comparison with GBHs does 
not extend to other states? Such a conclusion would have 
a fundamental impact on our understanding of the accre- 
tion process onto black holes and so it is very important 
to confirm, from timing properties, a definite non-soft state 
AGN. 

The only remaining case where there is at least some 
evidence for a non-soft state PSD is Ark 564. From ob- 
servations with RXTE, evidenc e for a low fre q uency break 
in the PSD was presented by [Pounds et al] (|200ll ). Evi- 
de nce for a second b r eak, a t high frequencies, was presented 
by iPapadakis et al. (I2OO2I) from analysis of ASCA obser- 
vations. iPapadakis et al.'n2002f ) interpreted the combined 
RXTE and ASCA PSD as a hard state PSD but did note 
that the high frequency break does not scale linearly with 
mass to the hard state of the archtypal GBH Cyg X-1; the 
frequency difference is too small. However as Ark 564 is 
one of the highest a ccretion rate AGN known (m_B ~ 1, 
iRomano et al.l (|2004l ')'). a VHS state remains a strong pos- 
sibility. 

Unfortunately proper modelling of the combined long 
and short time-scale PSD is hampered by the gaps which 
occur in the ASCA light curves, and hence in the PSD, at 



the orbital period (~ 5600s). We therefore made a 100 ks 
continuous observation with XMM-Ncwton to fill that gap 
in the PSD. Spectral analysis of this XMM-Ncwton obser- 
vation has been present ed by IPapadakis et al.l (|2007h and 
iBrinkmann et alT (|2007l ) and cross-spectral analysis (i.e. 
lags between energy bands as a functi on of Fourier fr e quenc y 
or time-scale) has been presented bv I Arevalo et al.l (|2006l ). 

Here we present a full power spectral analysis of the 
combined RXTE, ASCA and XMM-Newton data and find 
firm evidence for two breaks in the combined PSD. We also 
show that the PSD can be parameterised as the sum of two 
Lorentzian components. We note particularly that the high- 
est frequency Lorentzian, for which good spectral informa- 
tion is available from both the XMM-Newton and ASCA ob- 
servations, is stronger at higher energies, exactly as is found 
for Lorentzian components in GBHs. In GBHs in the hard 
and VHS states, changes in the lag between bands occur 
at approximately the time-scales at which the PSD changes 
from one Lorentzian component to another, supporting the 
model that each Lorentzian comes from a separate physical 
emission region. In this paper we therefore re-examine the 
lag spectrum and compare it with the PSD, to test whether 
a two-component emission region may apply in Ark 564 and 
also as a further diagnostic of the state of the AGN. 

The paper is arranged as follows: In Sec. [2] we describe 
the observational data sets which we used and discuss the 
individual ASCA and XMM-Newton PSDs in Sees |3] and H 
We then discuss the joint RXTE ASCA and XMM-Newton 
power spectrum and its energy dependence in Sec. [S] and 
compare the powerspectrum with the lag spectrum in Sec|6l 
We summarise our conclusions in Sec. [T] 



2 SUMMARY OF OBSERVATIONS 

2.1 RXTE 

Ark 564 was observed by RXTE from 1999 January 1 to 
2003 March 4. For most of this period observations were 
carried out approximately once every 4 days but from 2000 
May 31 to 2000 July 01, observations were carried out 8 
times per day. The observations from 1999 January 1 to 
2000 S eptember 19 are described in detail in [Pounds et al.l 
(|200ll 'l. We have extracted all of the RXTE data from the 
archive and produced 2-10 keV lightcurves in the same man- 
ner as we have descri bed for other _RXTE AGN monitoring 
programmes (e.g. see lM'^Hardv et aLll2004l . for details). For 
use in power spectral analysis we produce a long time-scale 
light curve, covering the full observing period with 4-day 
sampling, and an intermediate time-scale lightcurve, cover- 
ing the one month period of 3-hour sampling. 

2.2 ASCA 

Ark 564 was monitored with ASCA between 2000 June 1 
and July 5, producing a month-long light curve. This obser- 
vation was only interrupted by the regular Earth occulta- 
tions of the satellite on the orbital time-scale of 5632s. Bin- 
ning the data in orbit-long bins produces an evenly sampled 
light curve, containing 551 points, that probes time-scales 
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from a few times 10 to 10 s. Additionally, each of the 
551 individual orbits cover time-scales of 10^ — 10^ s, whose 
power spectra can be combined to produce a high quality 
high frequency PSD. 

We downloaded SIS data from the TARTARUS 
database in 0.6-2 and 2-10 keV energy bands and pro- 
duced background subtracted light curves, with average 
count rates of 2.7 and 0.6 c/s respectively. The average 
photon energies of these two bands are 1.18 and 3.2 keV 
respectively. 

The lightcurves, PSD and other vari ability properties 
of this data s et have been studied by e.g. iPapadakis et al.l 
(,2007. '); ,Edelson et all (|2002l ). 



2.3 XMM-Newton 

Ark 564 was observed by XMM-Newton from 2005 January 
5, 19:47, to 2005 January 6, 23:16. The European Photon 
Imaging Cameras (EPIC) PN, MOSl and M0S2 cameras 
were all operated i n small window mode us ing a medium fil- 
ter, as described bv lPapadakis et all l|2007h . The total expo- 
sures for the PN and MOS cameras were 98.8 ks and 99.1 
ks respectively. For the PN camera, we extracted source 
photons from a ~ 2' x 2' region and the background was 
selected from a source-free region of equal area on the same 
chip. Only single and double events, with quality flag=0 
were retained. For the MOSl and M0S2 cameras, source 
photons were extracted from a circular region of ~ 46" in 
radius and single, double, triple and quadruple events were 
used. The PN data was free of pile-up, but both MOS cam- 
eras suffered significant pile-up in the cores of the target 
PSFs. Therefore we discarded data from the central 12" di- 
ameter of each MOS exposure. The average PN background- 
subtracted count rate in the 0.2-10 keV band was ~ 28 c/s. 
The background level was low and almost constant through- 
out the observation. 

We produced light curves in the 0.2-2 and 2-10 keV 
energy bands in 48 s bins, and combined the light curves 
from the different detectors to produce weighted average 
light cur ves. The 0.2-10 keV l ightcurve is shown in Fig [T] 
(see also 'Papadaki s et aLlbOOTj ). As the variability proper- 
ties normally depend on energy band, we also constructed 
XMM-Newton light curves from the PN data alone in the 
0.5-2 and 2-8.8 keV bands, to match the average energies 
of the ASCA light curves. We used these lightcurves in all 
quantitative analyses. 



3 ASCA POWER SPECTRUM 

The superposed PSDs derived from the ASCA 2-10 keV and 
0.6-2 keV data are shown in Fig. [S] Data below 10~''Hz re- 
sult from the light curves binned on the orbital time-scale. 
Data above 10~^'^IIz are the sum of the PSDs derived from 
the many short sub-orbital light curves. Note the frequency 
gap that these data do not cover. The Poisson noise level 
has been subtracted. On time-scales longer than the orbital 
time-scale, the low and high energy PSDs have a similar 
normalisation. However on shorter time-scales the 2-10 keV 




Figure 1. A weighted average combined PN, MOSl and MOS2 
XMM-Newton lightcurve in the 0.2-10 keV band. 




-4 -3.5 
Log Freq (Hz) 

Figure 2. ASCA PSDs. The open (black) squares, with thick 
lines, are the 2—10 keV data. The filled (red) circles, with 
thin lines, are the 0.6—2 keV data. The Poisson noise level has 
been removed from both datasets. Data below 10~^Hz result 
from lightcurves binned on the orbital time-scale. Data above 
are the sum of the PSDs derived from the many short 
sub-orbital lightcurves. The 2-10 and 0.6-2 keV data points have 
been shifted very slightly in frequency to avoid overlap at both 
ends of the overall PSDs. 



PSD lies above the 0.6-2 keV PSD, except perhaps at the 
very highest frequencies (~ lO'^^Hz) where the two PSDs 
may be converging. This frequency dependent excess of the 
high energy PSD over the low energy PSD both below and 
above the high frequency break (~ 10~^Hz), which is dif- 
ferent from the behaviour of other well observed AGN, (e.g. 
NGC4051, M'Hardv et al.ii2004 ). will be modelled properly 
in Section [S] 



4 XMM-Newton POWER SPECTRUM 

The XMM-Newton light curves were almost completely 
continuous except for 5 small gaps of duration 20-30 sees 
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Figure 3. XMM-Newton PSDs. Here, for display purposes, we 
show the 2—10 keV (upper data points, black open squares, thick 
lines) and 0.2—2 keV (lower, red filled circles, thin lines) PSDs, 
with Poisson noise level subtracted, as the band centroids are 
more separated than for the 0.5-2 and 2—8.8 keV bands. Thus 
spectral differences, particularly the excess around 10~^Hz at 
higher energies, are a little more pronounced. We do, however, 
use the 0.5-2 and 2-8.8 keV bands (Table [TJ for quantitative 
analyses, to match better the ASCA energy bands. 




io-> 

Frequency [Hz] 



Figure 4. XMM-Newton 2-8.8 keV XMM-Newton PSD showing 
the model components. The broad high frequency Lorentzian is 
shown by a dash (light blue) line, the narrow Lorentzian is shown 
by a dot-dash (green) line, the underlying power law component 
is a dotted (black) line, the data points are (red) crosses, and 
the overall model is a solid (red) line. The Poisson noise level has 
been subtracted. 



each. Following the procedure established in lM'^Hardv et all 
|2004) these gaps were filled by linear interpolation, to- 
gether with the addition of noise characteristic of the sur- 
rounding data points, in order to eliminate spurious high 
frequency power in the power spectrum. Power spectral den- 
sities were derived from these light curves using a standard 
direct Fourier transform. In Fig. |3]we show the 0.2-2 and 
2-10 keV PSDs and we see close agreement with the ASCA 
PSDs in the same energy bands. As with the ASCA PSDs, 
the low and high energy XMM-Newton PSDs coincide very 
closely at frequencies below 10~*Hz, but at higher frequen- 
cies the 2-10 keV PSD exceeds the lower energy PSD. Note 
that here, for the low energy band, we show the 0.2-2 keV 
PSD rather than the 0.5-2 keV PSD and note that the bump 
centred around 10~^Hz is even less pronounced in the 0.2-2 
keV band, strengthening our observation that the excess has 
a strongly energy-dependent shape. Although not shown, 
the XMM-Newton 0.5-2 keV PSD is almost identical to the 
ASCA 0.6-2 keV PSD. 

A simple power law fit to both low and high energy 
PSDs is a very poor fit, leaving a large residual bump cen- 
tred around 10~"^Hz. A bending powerlaw, however, which 
provides a good fit to the high/soft state of Cyg X-1 and 
to th e XMM-Newton PS Ds of NGC4051 llM'^Hiirdv et al.l 
I200I) and MCG-6-30-15 jWRardv et"alll2005l ). is a good 
fit (Table [1} . Tying the break frequencies to be the same 
in both energy band PSDs, we note (as can also be seen 
directly from Fig. |3]), that although the slopes above the 
break are similar in both energy bands, the slope below the 
break is noticeably flatter at higher energies. The ASCA 
PSDs, although less well defined because of the gap on 
the orbital time-scales, are quite consistent with the fit pa- 
rameters given in Table [T] This behaviour is different from 
that of NGC4051 and MCG-6-30-15 (e.g. see Table 6 of 



iM'^Hardv etHlbOOSh where, although the slopes below the 
break are not very well determined, there is no evidence for 
a variation with energy. Above the break, however, the slope 
is flatter at higher energies in NGC4051 and MCG-6-30-15. 

The PSD slope variations seen here in Ark 564 can be 
more naturally interpreted in terms of a power law (which 
is necessary to explain the power at low frequencies) to- 
gether with a Lorentzian component which is stronger at 
higher energies. Such a model, in which the power law slope 
and normalisation, and Lorentzian central frequency are 
tied, provides an equally good fit to the data. In the case 
of the latter model, we note that the Lorentzian compo- 
nent is approximately three times stronger in the 2-8.8 keV 
than 0.5-2 keV band (Table [2|. An increasing Lorentzian 
amplitude at higher energies is also typically observed in 
QP Os and some broad Lor entzian components in GBHs 
(e.g. lPottschmidt et al.ll2003l ). As the exact values of the fit 
parameters depend on proper fitting of the lower frequency 
data from RXTE and ASCA we leave further discussion to 
Section [5] 

There is a residual to the power law plus Lorentzian fit 
to the XMM-Newton PSD at ~ 2.2 x lO'^Hz. This residual 
can be fitted (using the unbinned data) by the addition of a 
second, narrow (width ~ 8 x 10~''IIz), Lorentzian with nor- 
malisation approximately half that of the broad Lorentzian. 
The fit is shown in Fig. 3] The fit is improved but, adjusting 
the F-test probabilities to take account of the fact that the 
line could have appeared in any of the available frequency 
bins, we find that the improvement is significant only at the 
90% confidence level. 
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Table 1. Bending Powerlaw: combined fit to the XMM-Newton soft and liard band data 

Energy Band a^ia oihigh I'B red d.o.f 

keV Hz 



0.5-2 1.31 ±0.12 3.62 ±0.37 1.73 ± 0.44 x 10"^ 

2.0-8.8 0.95 ±0.14 3.57 ± 0.61 tied 0.88 37 



Table 2. Power Law with Lorentzian Line: combined fit to the XMM-Newton soft and hard band data 

Energy Band a Line centre Line width Line normalisation red d.o.f 

keV Hz 



0.5-2.0 1.80 ±0.19 3.94 ± 1.7 X 10-'' 6.8 ± 16 x lO"* 8.8 ± 5.7 x 10"^ 

2.0-8.8 tied tied 1.1 ± 0.17 x 10"^ 24.1 ± 8.4 x lO'^ 1.05 35 



5 COMBINED XMM-Newton, ASCA AND RXTE 
POWER SPECTRA 

The RXTE data, which covers time-scales fr om ~ 3 years 



to 3 hours, has aheady been analyzed by [Pounds et al 



2001 



and, including further data, by iMarkowitz et al' 



20031) who both fit a breaking power law model to the 



PSD. IMarkowitz et al.l (|2003l ) derive a slope at the low- 
est frequencies ol of 0.05^2 05' * break frequency, vl = 
1.591q'95 ^ 10~^Hz and an intermediate slo pe above the 
break of ai = 1.20t^;3^. The earlier results of lPounds et al.l 
l|200lh are consistent with these values. We begin by refit- 
ting these data using o ur standard simula tion-based mod- 
elling method (e.g. see lUttlev et al.l 



20021). We fit a sim- 



M'^Hardv et"alll2004 



pie b ending power law model (e.g. see[ 
I2OO5I ). We bin the long time-scale RXTE light curve into 4- 
d bins, and the short time-scale RXTE light curve into 3-h 
bins. The resultant best -fitting values are similar to those of 
IMarkowitz et all (|2003h and the errors in the fit parameters 
are similarly large. As our value of ol is quite consistent 
with zero and as, in both their hard and very high states, 
GBHs have ah = 0, we fix ol = in subsequent power law 
fits. 



5.1 Doubly-bending power law fit 

We then fit the combined RXTE, ASCA and XMM-Newton 
data sets with a doubly bending power law model, including 
a second, higher frequency break at uh and slope an at the 
highest frequencies. We include the ASCA 2-10 keV SIS 
data in 512s bins and the XMM-Newton PN data in 48s 
bins in the 2-8.8 keV band to match the average energy of 
the ASCA band. The PN bin value was chosen so that the 
resultant light curve has 2^' points so that we can use the 
fast Fourier transform, rather than the slower direct Fourier 
transform, for that part of the simulation, without loss of 
accuracy. 

We obtain a very good fit (acceptance probability 
P=75% - see Table E) and Fig. ^ with ai = 1.2toi, 
VL = T.Stg^s^ x lO-'^Hz and vh = 2.4^^:^ x 10"^ Hz. These 
break frequ encies are in good agr eement with the values 
obtained bv lPapadakis et al.l l|2002l ) from analysis of ASCA 
data. 

In order to make a combined PSD fit of the soft X-ray 
band, we first compared the low frequency variability power 
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Figure 5. PSD of Ark 564 from RXTE, ASCA and XMM- 
Newton data in the 2-10 keV RXTE energy band. The solid 
black line shows the best-fitting double-broken power law model. 
The symbols with error bars represent the data, unfolded from 
the distortions of the observational sampling pattern. The crosses 
(blue) are from the RXTE 4-d sampling lightcurve, the asterisks 
(green) are from the RXTE 3-h sampling lightcurve, the open 
squares (black) are from the ASCA 2-10 keV SIS 512s-binned 
lightcurve and the filled circles (red) are from the 48s-binned 
XMM-Ncwton 2-8.8 keV PN data. 



of the ASCA SIS data in the 0.6-2 and 2-10 keV bands. 
The resultant power spectra showed identical amplitudes, 
so both energy bands can be fit together at low frequen- 
cies without needing to rescale their PSDs. We assume that 
the close similarity of the 0.6-2 and 2-10 keV lightcurves 
continues to even lower frequencies, i.e. those sampled by 
the RXTE long term data. We therefore fitted the RXTE 
2-10 keV long term data together with the SIS and XMM- 
Newton 0.6-2 keV data to cover the same frequency range 
as above. The fit is again very good (P=85% - TableO giv- 



ing parameters aj = 1.5 



+0.1 



Uh = 1. 



j+l.l 



3.6^2.9 X 10" 



Hz and 



X 10 Hz. The steeper intermediate slope re- 



flects the lower high-frequency power in the soft band com- 
pared to the hard band. The break frequencies in both bands 
are consistent within 1 cr uncertainties, as can be seen from 
the confidence contours plotted in Fig. |6] 
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Figure 6. 68 % and 90 % confidence contours for high and low 
bend frequencies of the doubly-bending model fit to the combined 
PSD. Hard band contours (black and red respectively) are plotted 
in solid and dashed lines, soft band contours (green and blue 
respectively) in dot-dashed and dotted lines. 



5.2 Double Lorentzian fit to the power spectrum 

Motivated by the strong evidence for a Lorentzian compo- 
nent in the XMM-Newton data alone, we have fitted the 
combined RXTE, ASCA and XMM-Newton data with a 
model consisting of two Lorentzian components. 



-At 



<H+4QH(«^-i^c,H) 



where the Lorentzian peak frequency Vpc^k is related to the 
centroid frequency Vc by Vpeak = fc(l + l/iQ^y^^ and the 
subscripts L and H refer to the low and high frequency com- 
ponents, respectively. We fitted the combined PSD in soft 
and hard energy bands independently, allowing all param- 
eters to vary. The low frequency component converged to 
very low values of the quality factor Q = Vc/ A^^vvhm , i.e. to 
very wide Lorentzians where the Lorentizan shape is almost 
insensitive to the value of Q. We therefore fixed this param- 
eter at Q = 0.01 (with 90% and 99% upper confidence limits 
of Q = 0.2 and 0.5 respectively) for both energy bands to 
enable a more detailed parameter grid search for the other 
parameters, and hence to obtain better constraints on these 
parameters. The best-fitting values are listed in Table [3] 
The Lorentzian peak frequencies are similar to, but not ex- 
actly the same as, the corresponding break frequencies, with 
the Lorentzian frequencies lying closer to the centre of the 
band-limited power. 

Figures [7] and [S] show the best-fitting two-Lorentzian 
models to the hard and soft band respectively. In this inter- 
pretation, the increase in variability power at high frequen- 
cies in the hard band is entirely due to the greater strength 
of the high frequency Lorentzian in this band. The peak 
frequencies in both energy bands are consistent within the 
90% confidence limits (see contour plot in Fig[5}. 

Note that Figs.[7]and[S]show the PSD unfolded through 
the distortions produced by sampling effects. These distor- 
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Figure 7. PSD of Ark 564 from RXTE, ASCA and XMM- 
Newton data in the 2-10 keV RXTE energy band. The dashed 
(red) line shows the low frequency Lorentzian, the dotted (blue) 
line shows the high frequency Lorentzian and the solid black line 
shows the best-fitting combined double Lorentzian model. The 
symbols with error bars, which are the same as in Fig. |5] repre- 
sent the data unfolded from the distortions of the observational 
sampling pattern. 



tions also depend on the underlying PSD shape so the re- 
sultant unfolded PSD depends on the model being fitted. If 
the two-Lorentzian model is the correct underlying PSD, the 
strong power component around 10~^Hz will lead to spuri- 
ous apparent power at lower frequencies in a PSD which has 
r(l) not been unfolded from the distortions, and where the low 
frequency part of the PSD has been derived from observa- 
tions which do not sample the 10~^Hz region. This effect is 
known as aliasing (e.g. see lUttlev et al ] |200l for more de- 
tails) and here particularly affects the long timescale RXTE 
data (blue crosses in Figs [S] [7] and [Sj . When we unfold the 
PSD we remove this spurious low frequency power. Thus the 
lowest frequency part of the PSDs in Figs. [7] and |8] lies be- 
low the same part in the bending power law model (Fig[5]), 
where the model did not contain so much power at higher 
frequencies. Therefore the low break frequency in Fig. [5] is 
lower than in Figs. [7]and|8l 

Since the normalisations of low and high-frequency 
Lorentzians are highly correlated with each other and with 
the Lorentzian frequency, we only treat the ratio of the low 
and high-frequency Lorentzian normalisations, Al/Ah, as 
a free parameter and do not fit the high and low-frequency 
normalisations separately. Although the 1-dimensional er- 
rors quoted in Table 3 seem to indicate that the normali- 
sation ratios for the hard and soft band overlap, we show 
in Fig. [To] that the 90% confidence contours of i-'pcak,H ver- 
sus Al/Ah do not overlap. Therefore, for the same i-'peak,H, 
the high-frequency Lorentzian is always relatively stronger, 
compared to the low-frequency Lorentzian, in the hard band 
than in the soft band. In fact, the low- frequency Lorentzian 
normalisation is similar in both bands and relatively well- 
constrained. For the best fitting model parameters, the val- 
ues of fractional rms in the low frequency Lorentzian are 
35 ± 1% and 33 ± 1% in hard and soft bands respectively, 
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Figure 8. PSD of Ark 564 from RXTE in the 2-10 keV band 
and ASCA and XMM-Newton data in the 0.6-2 keV energy band. 
The solid line shows the best-fitting double Lorentzian model and 
symbols with errorbars, which are the same as in Fig. |5] repre- 
sent the data, unfolded from the distortions of the observational 
sampling pattern. 
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Figure 9. 68 % and 90 % confidence contours for high and low 
peak frequencies of the two-Lorentzian model fit to the combined 
PSD. Hard band contours are plotted in solid black and dashed 
red lines, soft band contours in dot-dashed green and dotted blue 
lines. 



and in the high frequency Lorentzian the values are 15 ± 1% 
and 10 ±1% in hard and soft bands. Thus the low frequency 
Lorentzian rms values are consistent with being the same, 
while the high frequency Lorenztian is weaker in the soft 
band, which explains the different hard and soft band PSD 
shapes at high frequencies. 



6 IMPLICATION OF LORENTZIAN PSD 
SHAPE ON TIME LAGS BETWEEN 
ENERGY BANDS 

In Galactic black hole systems, lags are observed between 
light curves in different energy bands, with the higher en- 
ergy variations generally lagging the lower energy variations 
(referred to here as a 'positive' lag). If the hght curves in the 
different bands are separated into different Fourier compo- 
nents, it is generally found that the time lags between bands 
vary as a function of the Fourier frequency, i.e. time-scale, 
of the component, wit h larger lags being associated with 
longer time-scales. (See lNowak et all l| 19991 ) for a full discus- 
sion of the relevant analysis techniques) . Such lags have now 
been measured in some AGN (e.g Papadaki s ct al. 200li; 
M'^Hardv et"aLll2004l l. including in Ark 564 (|Arevalo et al.l 
20061). 

lArevalo et al] l|2006l ) show that, for Ark 564, the pos- 
itive time lag between the 0.5-2 and 2-8.8 keV XMM- 
Newton bands is approximately 800 s at low Fourier fre- 
quencies and ~ Os at high frequencies, with a sharp drop 
in time lag values at a frequency of ~ 10~* Hz. This 
stepped appearance of the lag spectrum resembles that seen 
in Cyg X-1 in the low/hard and intermediate states, where 
the steps in the lag spectrum correspond to the transi- 
tion frequencies between peaks in the corresponding PSDs 
(|Nowak et al.lll999l : lNowakll2000l '). The correspondence be- 
tween lag spectrum steps and PSD peaks suggests that each 
constant-lag value is associated with a given variability com- 
ponent, seen as a Lorentzian component in the PSD. There- 
fore, the shape of the lag spectrum in Ark 564 suggests that 
the two-component interpretation of the PSD is correct. 
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Table 3. Fit to the PSD using combined XMM-Newton, ASCA and RXTE data 

I'L [Hz] uu [Hz] aL 



Doubly-bending power law 

Energy Band 



aj 



fit probability 



0.6-2.0 
2.0-8.8 



3. 6+29 ^ 10-^ 
7.5128-1 X 10-^ 



1, 



+ 1-8 

-1.4 



X 10-3 

0.9 X 



0.0 
0.0 



1.5 
1.2 



+0.1 
0.4 
+0.2 
0.1 



3.4; 

4.2+ 



0.85 
0.75 



Two-Lorentzian model 

Energy Band 



low-!/ Lorentzian 
I'pcak.L [Hz] Ql 



high-;^ Lorentzian 
!^pcak,H[Hz] Qh 



^l/^h fit probability 



0.6-2.0 
2.0-8.8 



g.stj-g X 10- 

-+2.0 



9.5"* 



X 10" 



0.01 
0.01 



3.9 
6.1 



+2.2 
-2.3 
-1.5 



X 10 
X 10 



0.5^ 



-+0.5 
-0.38 

0.125+°-3** 



17.f 



.15 



+290.4 
-11.9 
+4.7 
1.0 



0.37 
0.40 



* The errors are large and are not very well defined. 
All errors are 68% confidence (Icr). 



6.1 Fitting the two-Lorentzian model jointly to 
the PSD and lag spectra 

To test the above hypothesis we directly fitted the lag spec- 
trum with a model corresponding to the lags expected from 
two Lorentzian variability components, each with a single 
distinct time lag, which is constant for all Fourier frequen- 
cies. The observed lag at any particular frequency is then 
a function of the overlap between the two Lorentzians, and 
must be calculated in the complex Fourier domain, by evalu- 
ating the cross-spectrum as a function of Fourier frequency, 
V. The cross-spectrum is given by C{u) — S{v)H*{v), where 
H{v) and S{v) are the Fourier transforms of the hard and 
soft band light curves and the asterisk denotes the complex 
conjugate. S{v) = Si{v)^-Sh{v) where the subscripts I and h 
denote the low and high-frequency Lorentzians respectively, 
and similarly for the hard band, H{v) — Hi{v) +Hh{v)- For 
the low and high frequency Lorentzian components, the soft 
band light curve is correlated with the hard band light curve 
except with a phase lag 4>i{v), <f>h{v) respectively. We deter- 
mine the phase lag from the fixed low and high-frequency 
Lorentzian time lags r; and th, by using 0(1^) = 2tivt. To 
avoid clutter, we now drop the frequency- dependence of the 
various parameters and write the real part of the cross- 
spectrum as follows: 



Re{C) = ^6(5") X Re{H) + Im{S) x Im{H) 
= [Re{Si) + ReiSh)][Re{Hi) + Re{Hh)] 
+ [Im{Si) + Im{Sh)][ImiHi) + Im{Hh)] 
- n + + cos(<;!.i + (l>h)] 

M + \Hh\cos{(j)i)] 
h + <t>i) + \Sh\sm{(t)i + <t)h)] 



= [15,1 cos( 
X [|Hi|cos( 
+ [|5',|sin( 



X [|H,|sin((^o) + |Hh|sin((^i)], 



(2) 
(3) 

(4) 



where Re and Im denote real and imaginary components. 
00 and <j)\ axe simply random phase values for the low and 
high frequency Lorentzians respectively, measured for each 
realisation of the light curves, and will be different for dif- 
ferent realisations. They are not related to the phase lags 
of interest between bands, i.e. and (j)h for the low and 
high frequency Lorentzians. For a noise process, such as the 
broad Lorentzians observed here, <jio and <j)-i, in indepen- 



dent Fourier frequencies, are independently and identically 
distributed over a uniform distribution between and 27r. 
The coefficients of the cross terms and will 

average to when taking the average value of the cross 
spectrum for many realizations of the light curves, or many 
frequencies within a frequency bin, and the average value of 
Re{C) is then 

{Re{C)) = \Si\\Hi\{cos,<j}i) + \Sh\\Hh\{cos<f,h) (5) 

Similarly we can show that Im(C) is given by 

{/m(C)) = \Si\\Hi\{sm(j>i) + |5h| liy^] (sin <^h> (6) 

1 5; I, \Hi\, \Sh\ and \Hh\ correspond to the square root 
of the power of those Lorentzians in each energy band, at 
the frequency bin to be evaluated. The observed phase lag 
is simply the argument of the cross-spectrum, so the ob- 
served frequency-dependent time lag, r(i/)ot)s can be evalu- 
ated from equations [5] and [6] using: 



1 - 
T(u)obs = - — tan 
2tvu 



ImjCju)) 
Re{C{u)) ■ 



(7) 



which can be fitted directly to the observed lag spectrum 
using the best-fitting Lorentzian PSD model obtained in 
Section 5.2 and letting the low and high-frequency time lags 
be free parameters. This simple model reproduces well the 
shape of the observed lag spectrum, shown in Fig. 1111 with 
a reduced = 1-35 for 17 d.o.f (fit probability=0.15). The 
best-fitting values are 638 ± 54 s for the low-frequency lag 
and —11.0 ± 4.3 s for the high-frequency lag. 

The low frequency components fitted to the PSD data 
in both energy bands are identical, but the high frequency 
components had different Q and peak frequency values, al- 
though they agree at the la level. Forcing the high fre- 
quency components to be equal in both energy bands, us- 
ing Qh = 0.5 and i^pcak.h = 6.1 x 10~^ Hz (both within la 
of the best-fitting values), produces a better fit to the lag 
spectrum, with reduced = 1-13 for 17 d.o.f. (fit probabil- 
ity = 0.32), and lags consistent with the previous fit. The 
consistency between the lag spectrum and PSD significantly 
strengthens the interpretation of the PSD of Ark 564 as the 
summation of Lorentzian components, and hence as a very 
high state system. 
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Frequency (Hz) 

Figure 11. Top paneliLag of soft (0.5-2 keV) band to hard (2- 
8.8 keV) band, as a function of frequency, calculated using ASCA 
and XMM-JVewton data. The soUd line represents the lag spec- 
trum expected for a two-Lorentzian PSD model with the param- 
eters obtained by fitting the PSD, forcing the Lorentizan com- 
ponents of the hard and soft bands to be equal. Each Lorentzian 
is assumed to have a constant time lag, whose values were ob- 
tained by fitting the lag spectrum, resulting in 638 ± 54s (i.e. 
hard band lags) and —11.0 ± 4.3s (hard band leads) for the low 
and high frequency Lorentzians, respectively. Bottom panehhard 
band PSD as shown in Fig [3 with individual Lorentzian com- 
ponents overplotted. Note how the lags (top panel) are approxi- 
mately constant in the frequency range when any one Lorentzian 
component dominates, and how the lags rapidly change in the 
frequency range where the dominance of individual Lorentzian 
components crosses over. 



6.2 Coherence 

The variability in the soft and hard bands is relatively co- 
herent. Measuring the degree of linear correlation between 
energy bands as a function of Fourier frequency shows that 
the light curves remain coherent (coherence, 7^, ^ 0.8 where 
7^ = for uncorrelated hght curves and 7^ = 1 for per- 
fectly correlated hght c urves) up to frequenc ies of at least 
10"^ Hz (see figure 3 in lArevalo et al] (120061 ) ). The coher- 
ence is consistent with the two-component interpretation, 
even if the variability patterns of the low and high fre- 
quency Lorentzians are uncorrelated. For example, if one 
energy band is fully modulated by one of the Lorentzians, 
while at the same time-scales the other energy band re- 
sponds equally to both Lorentzians, the coherence would 
drop to a value of 7^ — 0.5. This extreme case never oc- 
curs in the present data, so the coherence can remain over 
50 % even if the Lorentzians are uncorrelated. The rela- 
tive contribution of the Lorentzians to each energy bands 
is most different at frequencies ^ 2 x 10"'' Hz, where both 
Lorentzians have similar power in the hard band (Figs[7]and 
IHI, but the high-frequency Lorentzian is 10 times stronger 
than the low frequency component in the soft band. In this 



case, the coherence is expecte d to drop to ~ 0.75, as mea- 
sured by I Arevaio et al.l (|2006l ) (their figure 3). 

It is possible that the low-frequency variability compo- 
nent that we fit to the PSD as a very broad Lorentzian ac- 
tually contains more components. The available data, how- 
ever, are not good enough to disentagle extra components 
in the low frequency PSD. Unfortunately the lag spectrum 
cannot help in this case as the long term monitoring was per- 
formed only by RXTE, which has a hard energy response. 
We are therefore unable to compute a lag spectrum between 
hard and soft bands as we have been able to on intermediate 
and short time-scales. 



7 DISCUSSION 

7.1 Summary of PSD and lag fits 

We have fitted the PSD of Ark 564 over the 10~® - 10"^ 
Hz frequency range, combining long term archival RXTE 
monitoring data, a month-long observation by ASCA cov- 
ering intermediate frequencies and a new 100 ks observation 
by XMM-Newton which covers the shortest time-scales. We 
confirm the existence of two breaks in the PSD, which make 
this the only AGN so far observed to have band-limited vari- 
ability. 

The PSD can be fitted either with a double-bending 
power law, where the variability power drops at low and 
high frequencies, or with a combination of two distinct vari- 
ability components, peaking close to the low and high fre- 
quency breaks, respectively. Using two Lorentzian-shaped 
components, we obtain a good fit in both the soft 0.5-2 
and hard 2-10 keV bands, with self-consistent Lorentzian 
parameters between the energy bands. We note that, when 
fitting the upper break in the PSD, the frequency of the 
fitted Lorentzian is a factor of 4 lower than the fitted bend 
frequency as, in any band-limited PSD, the Lorentzian fits 
will tend to move towards the centre of the band. This dif- 
ference should be considered when combining values of PSD 
break timescales derived from different fitting methods. In 
the Lorentzian interpretation of the PSD, the increase in 
variability power at high frequencies in the hard band is ex- 
plained by a strengthening of the high frequency Lor entzian, 
as IS seen m GBHs (e.g. iPottschmidt et al.ll2003l ). while 
the low frequency Lorentzian remains approximately at the 
same level. 

The two-component interpretation is strongly sup- 
ported by the lag spectrum calculated between these en- 
ergy bands using ASCA and XMM-Newton data (Arevaio 
et al. 2006) . The lag spectrum has a step-like shape, where 
the time lags are approximately constant at low frequencies 
and drop sharply to a different constant value at around 
10~* Hz. The drop in the lag spectrum occurs exactly at 
the frequency expected if each Lorentzian component has a 
single time lag value associated with it. This characteristic 
lag spectral shape therefore provides strong evidence that 
the different Lorentzian components in the PSD arise from 
distinct physical regions, e.g. different radii in the accretion 
disc. 
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7.2 Scaling Relationships 

The lag spectrum of Ark 564 has a very similar shape to 
that of Cyg X-1 in its intermediate state (see figure 8 of 
lArevalo et al.|[2006l ) , although the frequenci es are, of course, 
very different. However our recent 3D fit (|M'^Hardv et al.l 
120061 ) shows that PSD break frequencies scale not only in- 
versely with black hole mass, M, but also approximately 
with rriE- If we apply the same scaling to the lag frequen- 
cies, then Ark 564 and Cyg X-1 in its intermediate state 
coincide almost exactly. 

Ark 564 was included in o ur 3D fit, using the value of 
Tb given bv iPapadakis et all (|2002l ). and lay close to the 
best-fitting plane. The better measured value of Tb which 
we derive here, on the assumption o f a bending power law, is 
identical to the value of Tb given bv IPapadakis et all (|2002l ) 
thus Ark 564 remains consistent with our previous scaling 
relationship. 

All previous studies of the scaling between AGN and 
GBHs have concentrated on the high frequency part of the 
PSDs as the (scaled) low frequency part of most AGN PSDs 
is not well measured. On the basis of just one AGN where we 
are able to fit a low frequency Lorentzian component we are 
therefore extremely cautious in any comparison, and we also 
note that intermediate state GBH PSDs show many and 
varied shapes where the separation betwee n high and low 
frequency Lorentzians is not c onstant (e.g. lAxelsson et al.l 
[2005, 200^: iBelloni et al.1 [20051 ). However the separation of 
approximately 1.5 decades in frequency between the two 
Lorentzians seen here is not wildly different from typical 
separations seen in intermediate state GBHs. Thus some 
sort of mass (and probably also accretion rate) related scal- 
ing probably applies to the lower, as well as to the upper, 
PSD frequencies. 



7.3 Radio luminosity and 'state' 

The PSD and lag spectrum indicate a VHS/intermediate 
state, rather than a soft state, with a hard state being ruled 
out because of the very high accretion rate. Is this state clas- 
sification consistent with other state measures, e.g. radio lu- 
minosity? VHS/intermediate state GBHs can be quite lumi- 
nous, although tra,nsient , radio sources (e.g. GRS1915-I-105, 
iMiller- Jones et al.ll2005l ) whilst soft state GBHs are unde- 
tected in deep radio observations. 

Ark 564 (black hole mass M ~ 2.6 x lO^^M©, 
iBotte e^l. 2004) is considerably more radio luminous than 
AGN with well defined soft-state PSDs an d similar black 
hole masses, e.g. NGC4051 (1.9 x IO'^Mm. JPeterson et al 
2004|) and MCG-6-30-15 (4.5 x lO^'Mo, iM'^Hardv et"al 
2005). Both NGC 4051 (M^Hardy et al. in pr eparation: 



see also .Christopoulou et al.l [l997l ) and Ark 564 l|Lal et al.l 
l2004l : ISchmi~ et al.ll200ll ) contain flat spectrum VLBI cores 
for which Luad-Arksm ~ 600 x LRad-NOCiosi (assuming 
redshift-derived distances). These observations are certainly 
consistent with a VHS/intermediate state interpretation for 
Ark 564 although the present lack of radio detection of soft 
state GBHs prevents a more quantitative comparison. 

If the same disc-jet coupling applies in NGC4051 and 



Ark 564 and the radio luminosi ty of the jet scales as m^ '^ 
(absolute, not Eddington, units: [Blandford fc Koniglll 19791 ), 

then we would ex pect LB.ad-Ark564 100 X Luad-NGCWS,! 

(deriving m from luttlev fc M'^Hardvll2005h . Thus it is not 
yet clear whether the higher radio luminosity of Ark 564 
relative to NGC4051 is purely due to a higher accretion 
rate or to some difference in jet structure. 



7.4 Physical origin of the high frequency 
Lorentzian timescale 

If we take the Lorentzian peak frequency (~ 6 x 10~* Hz, 
i.e. 1666s, in the 2-8.8 keV band) as being the best measure- 
ment of the underlying physical process responsible for the 
sharp drop in power, we can consider what that timescale 
might correspond to. The dynamical time-scale around the 
Ark 564 black hole (2.6 x 10^ Mq) is considerably shorter 
than the Lorentzian peak time-scale, unless we assume that 
the time-scale arises far (few tens of Rg) from the black hole. 
There is evidence f rom the broad iron X-ray flourescence 
lines (jPabianl [ioosi ) that the accretion disc often reaches 
down to < 27?g, close to the innermost stable orbit of a 
maximally spinning black hole black hole. If the Lorentzian 
peak time-scale arises from the inner edge of such a spin- 
ning black hole disc, it may well correspond to the viscous 
time-scale (assuming a viscosity parameter, a = 0.1). In 
that case the ratio of disc scale height to radius, H/ R, = 
1.2, which is reasonable given the very high accretion rate. 



7.5 Combined interpretation of PSD and lag 
spectrum 

A possible interpretation of both the PSD and lag spec- 
trum is that the variability originates mainly in two dis- 
tinct regions. Note that by 'variability' we mean the source 
of the variations, not the source of the X-rays. The vari- 
ability may be related to accretion rate fluctuations in the 
accretion disc. It is reasonable to assu me that variabilit y 
time-scales are larger at larger radii (e.g. lLvubarskii|[l997^ . 
Once produced, variations would then propagate inwards 
and, on reaching the X-ray emitting region, would modu- 
late the emission. It is also reasonable to assume that the 
emitted energy spectrum becomes harder towards smaller 
radii so, for example, the bulk of the 0.5-2 keV emission 
will be emitted from a larger radius than the bulk of the 
2-10 keV emission. If the region responsible for most of the 
variability at low frequencies lies outside all of the X-ray 
emitting region, then the time lag between the soft and 
hard bands may represent the time taken for the fluctua- 
tions to propagate from the 'soft' to 'hard' emissio n region 
(e.g. see lKotov et al.ll200ll . [Arevalo fc Uttlcv"2006') . If aU of 
the frequencies in the low frequency variability component 
are produced in the same region, for example a narrow an- 
nulus of the disc, then all of those frequencies will have the 
same soft-hard lag. 

If the source of the higher frequency variations lies 
within the X-ray emission region, e.g. between the soft 
and hard emission regions, then the high frequency vari- 
ations will affect only the harder emission region and the 
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resultant X-ray variations will be stronger at higher en- 
ergies, as we observe. In this scenario the predominantly 
softer, outer emission region, which is subject only to the 
slower variations, may form part, or possibly all, of the 
so-called 'soft excess', i.e. the excess in the energy spec- 
trum below ~ 2keV above an extrapolation of a power law 
to low energies. This softer region, although lying outside 
the harder emitting region is, however, probably still quite 
close (ie within ~ IORq rather than ~ 100-Rg) to the black 
hole. This sc e nario is co nsistent with the o bserv ations of 
[Turner et all 1I2OOII) ari d iBrinkmann et all l|2007l ') of Ark 
564. Turner et al. show that, on ~ 1 day time-scales, 

the soft excess emission from Ar k 564 varies at l e ast as 
much as the power law emission. IBrinkmann et ahl |2003) 
show, by fitting a combination of bremsstrahlung (not nec- 
essarily a physical bremsstrahlung component but simply 
to fit the 'soft excess') and power law components to 500s 
sections of the XMM-Newton data discussed here, that the 
'bremsstrahlung' component at low energies (0.3-1 keV) is 
less variable than the power law component in the same 
band, and that the power law component at higher ener- 
gies (3-10 keV) is more variable still. However a full under- 
standing of how the soft excess is related to the rest of the 
X-ray emission, e.g. whether it might a ctually be caused 
by absorption l|Gierliiiski fc Dong |2006| ) requires detailed 
frequency-resolved spectrosc opy (e.g. see [Papadakis ot al. 
l2007l : lRevnivtsev et al.ll200ll ) for which the present data are 
insufficient. 

We finally comment on the negative phase lag, — 11.0± 
4.3s, at the highest observable frequencies. As the high 
frequency X-ray variations at lower energies will also be 
produced at small radii, e.g. as the low energy tail of the 
harder emission region, there should be very little lag be- 
tween the soft and hard emission, as we see. However the 
lag is negative, i.e. the hard band leads, in the higher fre- 
quency component. This observation suggests that the high- 
frequency variability in the soft band might arise from re- 
flection/reprocessing within ~ IRg of the illuminating hard 
X-ray source, implying the existence of cold gas close to the 
central X-ray source. 
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